
BELLE 



Belle Preprint 2012-15 
KEK Preprint 2012-7 



CO 

o 



X 



CN 



Measurement of Branching Fraction and First Evidence of 
CP Violation in B° -> af (1260)tt^ Decays 



J. Dalseno, 24 ' 43 I. Adachi, 8 H. Aihara, 47 D. M. Asner, 36 V. Aulchenko, 2 T. Aushev, 14 
A. M. Bakich, 41 A. Bay, 21 K. Belous, 13 B. Bhuyan, 9 A. Bozek, 32 M. Bracko, 23 ' 15 
O. Brovchenko, 17 T. E. Browder, 7 V. Chekelian, 24 A. Chen, 29 P. Chen, 31 B. G. Cheon, 6 
K. Chilikin, 14 R. Chistov, 14 I.-S. Cho, 53 K. Cho, 18 Y. Choi, 40 Z. Dolezal, 3 Z. Drasal, 3 
S. Eidelman, 2 J. E. Fast, 36 V. Gaur, 42 N. Gabyshev, 2 A. Garmash, 2 Y. M. Goh, 6 
H. Hayashii, 28 Y. Horii, 27 Y. Hoshi, 45 W.-S. Hou, 31 Y. B. Hsiung, 31 H. J. Hyun, 20 

T. Iijima, 27 - 26 K. Inami, 26 A. Ishikawa, 46 R. Itoh, 8 M. Iwabuchi, 53 Y. Iwasaki, 8 
T. Iwashita, 28 T. Julius, 25 J. H. Kang, 53 C. Kiesling, 24 H. O. Kim, 20 J. B. Kim, 19 
Y. J. Kim, 18 K. Kinoshita, 4 B. R. Ko, 19 S. Koblitz, 24 P. Kodys, 3 S. Korpar, 23 ' 15 
P. Krizan, 22 ' 15 P. Krokovny, 2 B. Kronenbitter, 17 T. Kuhr, 17 T. Kumita, 49 Y.-J. Kwon, 53 
S.-H. Lee, 19 J. Li, 39 J. Libby, 10 C. Liu, 38 Z. Q. Liu, 11 R. Louvot, 21 J. MacNaughton, 8 
D. Matvienko, 2 S. McOnie, 41 K. Miyabayashi, 28 H. Miyata, 34 Y. Miyazaki, 26 
G. B. Mohanty, 42 D. Mohapatra, 36 A. Moll, 24 ' 43 N. Muramatsu, 37 M. Nakao, 8 
Z. Natkaniec, 32 E. Nedelkovska, 24 C. Ng, 47 S. Nishida, 8 K. Nishimura, 7 O. Nitoh, 50 
S. Ogawa, 44 T. Ohshima, 26 S. Okuno, 16 P. Pakhlov, 14 G. Pakhlova, 14 C. W. Park, 40 
H. K. Park, 20 T. K. Pedlar, 54 R. Pestotnik, 15 M. Petric, 15 L. E. Piilonen, 51 M. Prim, 17 

K. Prothmann, 24 ' 43 M. Ritter, 24 M. R5hrken, 17 H. Sahoo, 7 Y. Sakai, 8 T. Sanuki, 46 
O. Schneider, 21 C. Schwanda, 12 A. J. Schwartz, 4 K. Senyo, 52 O. Seon, 26 M. E. Sevior, 25 
M. Shapkin, 13 V. Shebalin, 2 C. P. Shen, 26 T.-A. Shibata, 48 J.-G. Shiu, 31 A. Sibidanov, 41 
F. Simon, 24 - 43 P. Smerkol, 15 Y.-S. Sohn, 53 E. Solovieva, 14 M. Staric, 15 M. Sumihama, 5 
T. Sumiyoshi, 49 G. Tatishvili, 36 Y. Teramoto, 35 K. Trabelsi, 8 M. Uchida, 48 S. Uehara, 8 
Y. Unno, 6 S. Uno, 8 P. Urquijo, 1 Y. Usov, 2 P. Vanhoefer, 24 G. Varner, 7 C. H. Wang, 30 
P. Wang, 11 M. Watanabe, 34 Y. Watanabe, 16 K. M. Williams, 51 E. Won, 19 B. D. Yabsley, 41 
Y. Yamashita, 33 Z. P. Zhang, 38 V. Zhilich, 2 V. Zhulanov, 2 and A. Zupanc 17 

(The Belle Collaboration) 

1 University of Bonn, Bonn 53115 
2 Budker Institute of Nuclear Physics SB RAS and 
Novosibirsk State University, Novosibirsk 630090 
3 Faculty of Mathematics and Physics, Charles University, 121 16 Prague 
^University of Cincinnati, Cincinnati, Ohio 45221 
5 Gifu University, Gifu 501-1193 
6 Hanyang University, Seoul 133-791 
7 University of Hawaii, Honolulu, Hawaii 96822 
8 High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 



Typeset by REVTeX 



1 



9 Indian Institute of Technology Guwahati, Assam 781039 
10 Indian Institute of Technology Madras, Chennai 600036 
11 Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 100049 
12 Institute of High Energy Physics, Vienna 1050 
13 Institute of High Energy Physics, Protvino 142281 
1J{ Institute for Theoretical and Experimental Physics, Moscow 117218 
15 J. Stefan Institute, 1000 Ljubljana 
16 Kanagawa University, Yokohama 221-8686 
17 Institut fur Experimentelle Kernphysik, 
Karlsruher Institut fur Technologie, Karlsruhe 76131 
18 Korea Institute of Science and Technology Information, Daejeon 305-806 
19 Korea University, Seoul 136-713 
20 Kyungpook National University, Daegu 702-701 
21 Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne 1015 
22 Faculty of Mathematics and Physics, 
University of Ljubljana, 1000 Ljubljana 
23 University of Maribor, 2000 Maribor 
2Jt Max- Planck- Institut fur Physik, 80805 Miinchen 
25 University of Melbourne, School of Physics, Victoria 3010 
26 Graduate School of Science, Nagoya University, Nagoya 464-8602 
27 Kobayashi-Maskawa Institute, Nagoya University, Nagoya 464-8602 
28 Nara Women's University, Nara 630-8506 
29 National Central University, Chung-li 32054 
30 National United University, Miao Li 36003 
31 Department of Physics, National Taiwan University, Taipei 10617 
32 H. Niewodniczanski Institute of Nuclear Physics, Krakow 31-342 
33 Nippon Dental University, Niigata 951-8580 
3Jt Niigata University, Niigata 950-2181 
35 0saka City University, Osaka 558-8585 
36 Pacific Northwest National Laboratory, Richland, Washington 99352 
37 Research Center for Electron Photon Science, 
Tohoku University, Sendai 980-8578 
38 University of Science and Technology of China, Hefei 230026 
39 Seoul National University, Seoul 151-742 
40 Sungkyunkwan University, Suwon 440-746 
41 School of Physics, University of Sydney, NSW 2006 
^ 2 Tata Institute of Fundamental Research, Mumbai 400005 
^Excellence Cluster Universe, Technische Universitat Miinchen, 85748 Garching 
^Toho University, Funabashi 274-8510 
Jr5 Tohoku Gakuin University, Tagajo 985-8537 
^ 6 Tohoku University, Sendai 980-8578 
^Department of Physics, University of Tokyo, Tokyo 113-0033 
^ 8 Tokyo Institute of Technology, Tokyo 152-8550 
^ 9 Tokyo Metropolitan University, Tokyo 192-0397 
50 Tokyo University of Agriculture and Technology, Tokyo 184-8588 
51 CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 



2 



Yamagata University, Yamagata 990-8560 

53 Yonsei University, Seoul 120-749 
5 ^ Luther College, Decorah, Iowa 52101 

Abstract 

We present a measurement of the branching fraction and time-dependent CP violation param- 
eters in B° — > o] t (1260)7r =F decays. The results are obtained from the final data sample contain- 
ing 772 x 10 6 BB pairs collected at the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. We obtain the product branching fraction 

B(B° of (1260)ttF) x B(af (1260) ir^ir^) = (11.1 ± 1.0 (stat) ± 1.4 (syst)) x 10~ 6 , 

and an upper limit on the product branching fraction for a possible decay with the same final state 

B(B° a±(1320)vrF) x 0(a± (1320) n^ir*) < 2.2 x 10~ 6 at 90% CL. 

In a time-dependent measurement to extract CP asymmetries, we obtain the CP violation param- 
eters 

A C p =-0.06 ± 0.05 (stat) ± 0.07 (syst), 
C CP =-0.01 ± 0.11 (stat) ± 0.09 (syst), 
S CP =-0.51 ± 0.14 (stat) ± 0.08 (syst), 

representing time and flavor integrated direct, flavor-dependent direct and mixing- induced CP 
violation, respectively. Simultaneously, we also extract the CP conserving parameters 

AC =+0.54 ± 0.11 (stat) ± 0.07 (syst), 
AS =-0.09 ± 0.14 (stat) ± 0.06 (syst), 

which, respectively, describe a rate difference and strong phase difference between the decay chan- 
nels where the af does not contain the spectator quark and those where it does. We find first 
evidence of mixing-induced CP violation in B° — > a] t (1260)7r =F decays with 3. la significance. The 
rate where the af does not contain the spectator quark from the B meson is found to dominate 
the rate where it does at the 4.1<r level. However, there is no evidence for either time- and flavor- 
integrated direct CP violation or flavor-dependent direct CP violation. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 
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FIG. 1: Leading-order tree (left) and penguin (right) diagrams for the decay B° — > a 1 vr^, where the 
parentheses in the figure indicate the two possible decays of the B°. The d quark may become part 
of the 7r~ or (q~. In the penguin diagram, the x in V xb refers to the flavor of the intermediate-state 
quark [x = u, c, t). 



I. INTRODUCTION 

CP violation in the standard model arises from a complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [H, 0]. Mixing- induced CP violation in 
the B sector has been clearly observed by the BaBar 0j and Belle [3] Collaborations in 
the b — > ccs induced decay B° — > J /if) Kg, while many other modes provide additional 
information on CP violating parameters. 

Decays that proceed dominantly through the b — > uud transition are sensitive to the 
interior angle of the unitarity triangle (f>2 (a) = arg(— V t dV* b ) / (V u dV* b ) . The BaBar and Belle 
Collaborations have reported time-dependent CP asymmetries in these modes that include 
decays such as B° tt+tt" 0,0, p% T @,@] and p + p~ @,[Io|. 

This paper describes the measurement of the branching fraction and time- dependent 
CP violation parameters of the b — > uud channel B° — > af n T , shown in Fig. [TJ The left 
diagram shows the dominant first-order or tree process while the right diagram shows the 
leading second-order loop or penguin process. This analysis can be used to test the QCD 
factorization framework which has been used to predict the branching fraction and CP 
asymmetries of this decay channel [ill 13|. Similar to p ± 7r =F , the state a^7r T is not a CP 



eigenstate; rather, it is a flavor non-specific state with four flavor-charge configurations that 
must be considered: B°(B°) — > a^n T (l4j ]. The combined information of the B flavor and 
a\ charge allows the determination of additional information compared to CP eigenstates. 
It allows us to separate the cases where the d quark from the B meson, which does not 
participate in the interaction (spectator), becomes part of the af or ir ± . 

The decay of the T(4S*) can produce a B°B° pair that must be coherent, of which one 
(B^ cc ) may be reconstructed in the afn 11 final state. This final state does not determine 
whether the B^ ec decayed as a -B or as a B°. The other B meson in the event (B^ ag ), 
however, can be reconstructed in a final state that determines its 6-flavor, and therefore the 
flavor of the B^ ec at the time of the B^ ag decay. The proper time interval between B^ cc 
and -£>Ta g ' which decay at time t Rec and t Ta _ g , respectively, is defined as At = t Rcc — ^Ta g - 
For the case of coherent B°B° pairs, the time-dependent decay rate in the quasi-two-body 
approximation when -B^ag possesses flavor q (B°: q = +1; B°: q = — 1) and the cq possesses 



4 



charge c (af: c = +1; a 1 : c = — 1), is given by [l 



e -\At\/r B o 

V{At,q,c) = {1 + cAcp)— O-f 
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8r B o 



(5cp + cA<S) sin Am d At - (C C p + cAC) cos Am d At 



Here, r^o is the -B lifetime and Amd is the mass difference between the two mass eigenstates 
of the neutral B meson. This assumes CPT invariance, no CP violation in the mixing, and 
that the decay rate difference between the two mass eigenstates is negligible. The parameter 
Acp measures time and flavor-integrated direct CP violation, 



A 



CP 



Y{B ->> a+7T") - Y(B -»■ a^n 
T(B a+vr-) + T(B -)• a^vr 



(2) 



The parameter Sep measures mixing-induced CP violation, and Ccp measures flavor- 
dependent direct CP violation. The quantity AC measures the rate asymmetry be- 
tween the flavor-charge configurations where the a± does not contain the spectator quark 
(r[Z?° — > afir~} + F[B° — > a~[7i + }), and where it does contain the spectator quark 
(T[B° — > ajf 7r + ] + T[B° — > a^TT^]), while AS is related to the strong phase difference 
between these two processes, 6. These two parameters are not sensitive to CP violation. 
Sensitivity to 02 comes from this relation between four of the measured parameters, 



S CP ± AS = y/l - (Ccp ± AC) 2 sin(2#f ± ± 5). 



(3) 



A feature common to this channel and the other b — > uud modes mentioned earlier is that 
an effective angle, (f)\ , is determined rather than 2 itself due to the possible presence of ad- 
ditional loop contributions. In the limit that only the dominant tree amplitude contributes, 
no flavor- dependent CP violation is expected and <p < f- ± = 02- Fortunately, this inconve- 
nience can be overcome with bounds on Ad> 2 = (j) 2 — (j)^ determined using either an isospin 
analysis [lH or SU(3) flavor symmetry 



From these parameters, (j) 2 can be determined up to a four-fold ambiguity [14 
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(4) 



These results can also be transformed into more physically intuitive parameters that char 



acterize direct CP violation in decays with particular topologies [16 



A, 
A^ 



■(Acp + C CP + A C p AC) 

1 + AC + AcpCcp 
-A C p + C CP + A CP AC) 
-1 + AC + AcpCcp 



(5) 



which describe CP violation involving diagrams where the a± contains and does not contain 
the spectator quark, respectively. 

The BaBar Collaboration has performed a branching fraction measurement of B° — > 
a^n T with 218 x 10 6 BB pairs [l7( and a time-dependent CP violation measurement with 
384 x 10 6 BB pairs [l8j. These results are collected in Table [H From a subsequent study of 
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TABLE I: Summary of physics parameters for B° 
tion 



mm 



a 1 ir^ obtained by the BaBar Collabora- 



Parameter 


Value 


B(B° -> af vr^) x B(af (1260) -> t^ttM 1 ) 


(16.6 ± 1.9 (stat) ± 1.5 (syst)) x 10" 6 


-4cP 


-0.07 ± 0.07 (stat) ± 0.02 (syst) 




-0.10 ± 0.15 (stat) ± 0.09 (syst) 


5cp 


+0.37 ± 0.21 (stat) ± 0.07 (syst) 


AC 


+0.26 ± 0.15 (stat) ± 0.07 (syst) 


AS 


-0.14 ± 0.21 (stat) ± 0.06 (syst) 



B — > J^i(1270)-7r and B — > i^i(1400)7r, the BaBar Collaboration has also obtained bounds 



on |A02| using SU(3) flavor symmetry [19 



Two separate measurements are described in this paper. In Sec. [TTJ we briefly describe 
the data set and Belle detector. The branching fraction measurement is described in Sec. IHIl 
There, we explain the selection criteria used to obtain signal candidates and suppress back- 
grounds followed by the fit method used to extract the signal component. After this, the 
results of the fit are presented along with a discussion of the systematic uncertainties. In 
Sec. IIV[ these same issues are described again for the time-dependent CP violation mea- 
surement followed by our conclusions in Sec. |Vj 



II. DATA SET AND BELLE DETECTOR 



This measurement of the branching fraction and time-dependent CP violation parame- 
ters in B° af 7r T decays is based on the final data sample containing 772 x 10 6 BB pairs 
collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) col- 
lider [2(|. At the T(45) resonance (y/s = 10.58 GeV), the Lorentz boost of the produced BB 
pairs was /?7 = 0.425 along the z direction, which is opposite the positron beam direction. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprising of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere j2l[. Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVD1) was used 
for the first sample of 152 x 10 6 BB pairs, while a 1.5 cm radius beampipe, a 4-layer silicon 
detector (SVD2) and a small-cell inner drift chamber were used to record the remaining 
620 x 10 6 BB pairs 0. We use a GEANT-based Monte Carlo (MC) simulation to model 



the response of the detector and determine its acceptance [23 
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III. BRANCHING FRACTION MEASUREMENT 



A. Event Selection 



We reconstruct B° — > afix T where af — > (n + Tt )ir ± . The af decay proceeds mainly 



through the p ^ and /o(600)7r ± intermediate states 24j. We assume that the p ^ interme- 
diate state gives the dominant contribution and treat the /b^OO)?!" 111 contribution separately 
in the systematic uncertainties. Thus, the signal MC events for establishing the selection 
criteria are generated as B° — > afir T decays where af — > p ^. The af mass and width are 



taken to be m ai = 1.23 GeV/c 2 and T ai = 0.40 GeV/c 2 (24 . 

Charged tracks are identified using a loose requirement of distance of closest approach 
to the interaction point (IP) along the beam direction, \dz\ < 4.0 cm, and in the transverse 
direction, dr < 0.4 cm. With information obtained from the CDC, ACC and TOF, particle 
identification (PID) is determined with the likelihood ratio + Cj). Here, £j (£j) is 

the likelihood that the particle is of type i (j). A requirement of J~-k/-k < 0.4 is placed on all 
charged pion candidates, which retains 91% of all pions from B° — > af ir T , but only 9% of 
kaons. To further suppress background from particle misidentification, vetoes are applied on 
particles consistent with the electron or proton hypotheses. Additional SVD requirements 



of two z hits and one r — (f) hit |25j are imposed on the charged tracks so that a good quality 



vertex of the reconstructed B candidate can be determined. 

The intermediate dipion state is reconstructed above the Kg region with an invariant mass 
0.52 GeV/c 2 < m(7r + 7r~) < 1.1 GeV/c 2 . This is combined with another pion and forms an 
af candidate if the invariant mass is in the window 0.85 GeV/c 2 < m 37r < 1.75 GeV/c 2 , 
below the charm threshold. Upon combination with another pion, a B candidate is formed. 

Reconstructed B candidates are described with two nearly uncorrleated kinematic vari- 
ables: the beam-energy-constrained mass Mb c = a/ (-^bcam) 2 — (Pb MS ) 2 ; an d the energy dif- 
ference AE = E% MS - Eg^, where -E^m is the beam energy and E% MS Qo§ MS ) is the 
energy (momentum) of the B meson, all evaluated in the center-of-mass system (CMS). The 
B candidates that satisfy M^ c > 5.27 GeV/c 2 and \AE\ < 0.1 GeV are selected for further 
analysis. 

To reduce combinatorial background in forming the af candidate, the cosine of the angle 
between the prompt pion from af — > p°7r ± and the B candidate in the af rest frame is 
required to be between —0.85 and +0.85. This distribution is roughly flat for signal while 
peaking at ±1 for combinatorial background. This selection retains 80% of signal events 
while rejecting 43% of the combinatorial background. 

The dominant background in the reconstruction of -B^ ec is from continuum (e + e~ — > qq 
where q = u,d, s, c) events. Since their topology tends to be jet-like, in contrast to the 
spherical BB decay, continuum events can be distinguished from BB events using event- 



shape variables, which we combine into a Fisher discriminant J-'bi/qq |26[. The BB training 
sample is taken from signal MC, while the qq training sample comes from data taken below 
the T(4S') resonance. The Fisher discriminant is then constructed from the following 8 
variables: 

• | cos #tb,to I j where the angle is between the -B^ cc thrust direction and the thrust of 
the tag side. The thrust is defined as the vector which maximizes the sum of the 
longitudinal momenta of the particles. For a BB event, the pair is essentially at rest 
in the CMS, so the thrust axis of -B^ ec is uncorrelated with the thrust axis of -B^ag- I n 
a qq event, the decay products lie in the two jets which are back-to-back, so the two 
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thrust axes tend to be collinear. This variable provides the strongest discrimination 
against continuum. 

• \ cos 8 b, z \, where the angle is between the B^ ec flight direction and the z direction. 
This is the second most important variable in discriminating against continuum. BB 
pairs are produced in a correlated state, for which this variable follows a sine squared 
distribution, whereas for qq events, the distribution is uniform. 

• |cos0tb,z|j where the angle is between the S^ ec thrust direction and the z direction. 
This distribution is strongly influenced by detector acceptance at large values. As 
thrust and decay axes are related, this variable displays similar tendencies to the 
previously described variable. 

• Y2p? MS i where the transverse CMS momentum sum runs over all particles on the 
tag side. The continuum distribution generally has a higher mean because its decay 
product multiplicity is lower compared to BB. 

• Lq 11 = ^2 cn p™ S , where the CMS momentum sum runs over either the charged tracks 
(c) or neutral clusters (n) on the tag side. These distributions exhibit similar tenden- 
cies to the previous variable. 

• L2' 11 = J2cnPcJn S cos2 ^tef s ,TB; where the CMS momentum sum runs over either the 
charged tracks (c) or neutral clusters (n) on the tag side. The angle is between the 
particle direction and the £>^ ec thrust direction. In addition to the factors explained 
for the previous variable, the mean of the continuum distribution increases all the more 
due to higher values of cos # p cms TB . For jet-like events like continuum, the momentum 
of a particle in is closely aligned with its thrust, which itself is strongly correlated 
with the B^ cc thrust as explained for the first variable. 

The distributions for each of these discriminating variables are shown for simulated (MC) 
signal and continuum (data) events in Fig. El Before training, a loose requirement of 
cos#tb,to < 0.9 is placed that retains 90% of the signal while rejecting 50% of the con- 
tinuum background. The Fisher discriminant is also required to satisfy — 3 < Fbh/qq < 2. 

The next largest background comes from charm (b —> c) and charmless (b — y u, d, s) 
decays of the B meson and is found to exhibit peaking structure in the signal region due to 
the reconstruction of particular channels with a four-track final state. Defining the decay 
chain B° — y clyki, a\ — > p7r 2 , p — > 7r 3 7r4, we apply vetoes to remove these peaking backgrounds 
as summarized in Table HI1 

The background coming from b — > uud channels with the same final state as signal is 
studied separately. We consider a possible contribution from B° — y afir T by constructing 
the helicity variable H^-, defined as the cosine of the angle between the normal to the plane 
of the di candidate and the flight direction of the bachelor pion from the B evaluated in the 
3tt rest frame. 

On average, 1.6 B candidates are reconstructed per signal event. Selecting the best B 
candidate having the nearest M\> c with respect to the nominal B meson mass, the correct 
B is chosen in 91% of cases where the event contains multiple candidates. The fraction of 
mis-reconstructed signal events in general is 17%. As this procedure introduces bias to the 
M hc distribution, this variable is excluded from the fit to extract the signal yield. 
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FIG. 2: (color online) Simulated (MC) and off-resonance data distributions for the quantities used 
to construct the Fisher discriminant J~bi/ q g, normalized to have the same area. The solid blue 
histograms show signal MC events while the dashed red histograms show continuum events from 
data taken below the T(4S) resonance. 

The full selection used to define the event sample for the CP asymmetry measurement is 
also applied for the branching fraction determination. Since the B^ ec and -B^ag mesons are 
approximately at rest in the T(4S) CMS, the difference in decay time between the two B 
mesons, At, can be approximately determined from the displacement in z between the final 
state decay vertices, 

A* ~ fog ~ ^ Tag) = *L (6) 
Pic f3~fc K J 

The vertex of reconstructed B candidates is determined from the charged daughters 
using the known IP. The IP profile is smeared in the plane perpendicular to z to account 
for the finite flight length of the B meson in this plane. To obtain the At distribution, 
we reconstruct the tag side vertex from the tracks not used to reconstruct B^ jlBJ. The 
events must satisfy the requirements |At| < 70 ps and /iR CCj Tag < 500, where h Rcc Tag is the 
multi-track vertex goodness-of-flt, calculated in three-dimensional space without using the 
interaction-region profile's constraint. To reduce the necessity of also modelling the event- 
dependent observables describing the variable At resolution in the fit, the vertex uncertainty 
is required to be cr^ ec ' Tag < 200 /xm for multi-track vertices and cr^ ec ' Tag < 500 //m for single- 
track vertices. 
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TABLE II: Summary of peaking background vetoes. Alternate mass hypotheses have been applied 
to specific tracks where indicated in order to remove certain channels. The efficiency loss caused 
by these vetoes is also included. The X represents any charged track(s) that lead to a 4-body final 
state. 



Regions vetoed 


Modes vetoed 


Efficiency loss 


f .85 GeV/c 2 < m(7r 2 A 3 7r 4 ) < 1.89 GeV/c 2 
1.85 GeV/c 2 < m(7r 2 7r 3 A 4 ) < 1-89 GeV/c 2 


B -> D+fA'-vr+vr+lA 


3.7% 


3.06 GeV/c 2 < m(fj,xfi 2 ) < 3.14 GeV/c 2 






3.06 GeV/c 2 < m(/^i/x 3 ) < 3.14 GeV/c 2 




6.1% 


3.06 GeV/c 2 < m(nim) < 3.14 GeV/c 2 






0.480 GeV/c 2 < m(7r 2 7r 3 ) < 0.516 GeV/c 2 
0.480 GeV/c 2 < m(7r 2 7r 4 ) < 0.516 GeV/c 2 


B -> A^tt+tt^X 


11.9% 



After these selection criteria, the MC detection efficiencies are found to be 

SVD1 : e(B° -»■ af tt t ) = 0.1713 ± 0.0004 

SVD2 : e(B° -> af?r T ) = 0.2037 ± 0.0005, (7) 

where the uncertainty comes from limited MC statistics. Using independent data samples, 
we also determine correction factors to these efficiencies that account for the difference 
between data and MC. Correction factors in our reconstruction algorithm arise only from 
differences in P1D and are determined from an inclusive D* + — > D°[K~7r + }ir + sample to be 

SVD1 : rj(B° -> af tt t ) = 0.860 ± 0.031 

SVD2 : ri(B° -> af tt t ) = 0.855 ± 0.047. (8) 



We employ the flavor tagging routine described in Ref. [27J. The tagging information is 
represented by two parameters, the -B° ag flavor q and the purity r. The parameter r is an 
event-by-event MC determined flavor-tagging parameter that ranges from r = for no flavor 
discrimination to r = 1 for unambiguous flavor assignment. Due to a finite mistag probability 
w, the CP asymmetry is diluted by a factor 1 — 2w. The measure of the performance of 
the flavor tagging algorithm is the total effective tagging efficiency e c g = 6Ta g (l ~ 2u;) 2 , as 
the statistical significance of the CP parameters is proportional to (1 — 2u;) v /eTag, where 
exag is the raw tagging efficiency. These are determined to be e e s = 0.284 ± 0.010 and 
e cff = 0.301 ± 0.004 for SVD1 and SVD2 data, respectively. 



B. Event Model 

The branching fraction is extracted from a four- dimensional extended unbinned maxi- 
mum likelihood fit to AE, T^i^ m 37r and "H 37r from a data sample divided into 7 bins 
(/ = 0..6) in the flavor-tag quality r and 2 SVD configurations s. We consider 12 cate- 
gories in the event model: correctly reconstructed signal (referred to as truth signal here- 
after), mis-reconstructed signal, continuum, charm neutral and charged BB decays, charm- 
less neutral and charged BB decays, and five peaking backgrounds. In most categories, 
the linear correlations between fit variables are small, so the probability density function 
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(PDF) for each category j, is taken as the product of individual PDFs for each variable 
V*' S (AE, T b i /qq ,m 37n H 3n ) = V l ' s (AE) x V l > s (T bb/qq ) x V l > s {m 3n ) x V l > s (H 3n ) in each I, s bin, 
unless otherwise mentioned. 

The truth model shape is determined from correctly reconstructed signal MC events. 
The PDF for AE is taken to be the sum of two asymmetric-width (bifurcated) Gaussians 
incorporating calibration factors that correct for the difference between data and MC. These 
factors calibrate the mean and width of the core bifurcated Gaussian and are determined 
from a large-statistics control sample B° — > D~[K + 7t~7t~]tt + . The PDF for J^bl/qq here 
and throughout this analysis, for all 12 categories, is the sum of two bifurcated Gaussians 
in each flavor-tag bin /. The shape for all BB categories is fixed from the truth model 
except for the mean of the core distribution, incorporating calibration factors that correct 
for the shape difference between data and MC. The m 3jr distribution is modelled with an 
efficiency-corrected relativistic Breit-Wigner 

P„<™») = e (m3 , )j ^—^ J — m — ) , (9) 

where e s is the mass- dependent detection efficiency and T is the mass-dependent width 

Due to the finite width of the p°, the phase space p 1+s of the af decay, where the superscript 
represents the decay into a spin-1 meson in an 5*- wave configuration j28|, cannot be calcu- 
lated simply. Therefore, we model the phase space empirically with a 6th order Chebyshev 
polynomial 

6 

P 1+S (m 37r ) = 1 +^CjaK), (11) 
i=i 

where Ci is a Chebyshev polynomial of order i and q is the fit coefficient. The PDF for 7-L 3t: 
is a sum of symmetric Chebyshev polynomials up to 6th order. 

The mis-reconstructed model shape is determined from incorrectly reconstructed signal 
MC events. For AE, the PDF is taken as a smoothed one-dimensional histogram, the m 3lT 
PDF is the sum of two asymmetric-width Gaussians and the PDF for "% 37r is the sum of 
symmetric Chebyshev polynomials up to 8th order. 

The parameterization of the continuum model is chosen based on data taken below the 
T(4S) resonance; however, all the shape parameters are free in the fit to extract the branch- 
ing fraction. Since continuum is by far the dominant component, extra care must be taken 
to ensure this background shape is understood, so correlations above 2% are considered. 
The PDF for AE is taken to be a 1st order Chebyshev polynomial in each flavor-tag bin 
with a coefficient depending linearly on J 7 ^/^ 

V%(AE\J? a/ a) = 1 + (pJT + p\JF bl/qq )AE. (12) 

The m 3n shape is observed to shift quadratically with T h ij q ^ so the PDF is a sum of Cheby- 
shev polynomial up to 4th order that incorporates an offset. In addition, a small excess is 
seen above this distribution and is modelled with a Gaussian. The result is 



V l q l{m 3% \F b i /qq ) = (l-f s 



1 + Yl C kCk(m 3n -pt^bb/qq-P^bb/qq) 
k = l 



+ f s G(m 3n ;p s ,a s ). (13) 
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TABLE III: List of peaking backgrounds, assumed branching fractions and their respective expected 
yields after all selection criteria has been applied. 

Mode B (xlCT 6 ) SVD1 expected events SVD2 expected events 



B° -> p°p° 0.73 ± 0.28 2 11 

B° ->■ bflx^^]^ 0.16 ±0.03 5 23 

B° -> p°it + it- 4.35 ± 4.35 28 137 

B° -> vr+vr-vr+vr" 9.65 ± 9.65 26 125 

-> /op < 0.3 < 2 < 7 

B°^f fo <0.1 <1 <1 

£° -> /o7r + vr _ < 3.8 < 14 < 65 



The PDF for "H 37r is the sum of symmetric Chebyshev polynomials up to 6th order and a 
mirrored Gaussian 

+ f s G{U^ /A a s ) + f s G{U^ a s ). (14) 

The charm BB background shape is determined from a large sample of MC containing 
generic b — > c transitions and is subdivided into neutral and charged B samples. For 
AE, the PDF is a smoothed one- dimensional histogram; for 7713^, the PDF is the sum 
of Chebyshev polynomials up to 4th order; for H^, the PDF is the sum of symmetric 
Chebyshev polynomials up to 6th order plus a mirrored Gaussian. 

The charmless BB background shape is determined from a large sample of MC containing 
generic b — > u,d,s transitions and is subdivided into neutral and charged B samples. A 
sizeable correlation is seen between AE and and is taken into account with a smoothed 
two-dimensional histogram. The "H 37r PDF is the sum of symmetric Chebyshev polynomials 
up to 6th order with a mirrored Gaussian. 

Many decay channels contain the same final state as B° — > afix T . In addition to the 
possibility of B° — > afn T , we also consider those listed in Table ULTl which includes their 
expected yields in 772 x 10 6 BB pairs. With the exception of B° — > a^7r T , the assumed 
peaking background branching fractions are fixed in the nominal fit from the Heavy Flavor 
Averaging Group [29[. Where a mode is to be included in the fit model but only an upper 
limit is known, the branching fraction is taken as half the upper limit unless the mode 
contains an f , in which case it is assigned zero branching fraction and instead considered 
solely in the systematic uncertainties. The peaking background shapes are determined from 
individually generated MC samples and are fixed in the fit to data. The PDF for AE 
borrows the shape of correctly reconstructed signal events and includes a 1st order Chebyshev 
polynomial to model the mis-reconstructed contribution underneath. The m 3n PDF depends 
on the peaking background. For B° — > afir T , a sum of three Gaussians is used; for B° — > p°p° 
and bf n T , a smoothed histogram is used; otherwise, the sum of Chebyshev polynomials up 
to 4th order is used. The "H 3?r PDF is the sum of symmetric Chebyshev polynomials up 
to 8th order, except for B° — > p°p°, which is modelled with a smoothed symmetrized one- 
dimensional histogram. 



V s 



i + ]T 40,(^3.) 
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TABLE IV: Summary of yields fixed relative to other yields free in the fit where the uncertainties 
are from limited MC statistics. The mis-reconstructed yield is fixed relative to the signal yield, 
and the charm and charmless B + B~ background yields are fixed relative to their respective B°B° 
background yields. 



Yield 


SVD1 


SVD2 


JV Mis 
Archarm; s 

B+B~ 
» rcharmlcss; s 

B+B~ 


(0.167 ± 0.001) N§y m 
(1.585 ±0.019)A^ a ™ ;SVD1 
(0.568 ±0.006)^ h o a ™ lcss;SVD1 


(0.166 ±0.001)A^ D2 
(1.700 ± 0.010) A^ a ™ ;SVD2 
(0.574 ±0.003)iY^, a ™ lcss;SVD2 



The total likelihood for 208238 B° — > afn T candidates in the fit region is 

- £ - Ns E; f l ' s Nl >" 

£ = n e v, iie ^fM( AEi > ^ (is) 

l,s *> s ' i=l j 

which iterates over z events, j categories, / flavor-tag bins and s detector configurations. In- 
stead of two free signal yields , the branching fraction is chosen as a single free parameter 
and is incorporated into the fit with 

N° lg = B(B° -+ af^) x B(af nV^N^rj^, (16) 

where eg ig and 7/g ig are given in Eqs.[7]and[8], respectively. Similar conversions are done for the 
remaining peaking backgrounds using their expected values from Table IHIl The fraction of 
events in each flavor-tag bin I, for category j, is denoted by f^ s - The fraction of signal events 
in each l,s bin, / S j*, has been calibrated with the B° — > D~[K + tt~7t~]tt + control sample. 
Other free physics parameters include the af width and the product branching fraction 
B(B° -> afir*) x B{af -)• irVir*). A1so free are the yields N ^ and Nf^ less ' s ; 

the remaining yields are fixed to the values given in Table [IV] as determined from MC. In 
total, there are 121 free parameters in the fit. In addition, all shape parameters of the 
continuum model are free in the fit to data. 

C. Fit Result 

We perform a fit to the data, with the projections shown in Fig. [31 and obtain the product 
branching fractions 

B(B° -> a^(1260)vr T ) x B(af (1260) -> vr ± 7r T vr ± ) = (11.1 ± 1.0 (stat) ± 1.4 (syst)) x 10~ 6 , 

(17) 

corresponding to a yield of 1445 ±216 events, and 

B(B° -> af (1320)7T T )xi3(af (1320) ->■ 7r ± 7r T vr ± ) = (1.5±0.4 (stat)±0.4 (syst)) x 10^ 6 , (18) 

corresponding to a yield of 282 ± 106 events. The statistical correlation coefficient between 
these two measurements is —0.41 and the statistical significance of the af peak is 16a, 
estimated by comparing the likelihood of the nominal fit result with that of a fit where 
the -B — > afir T branching fraction is fixed to zero. We also measure the af width to be 
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FIG. 3: (color online) Projections of the fit to the B° — > a^it^ data. The points with error 
bars represent the data and the solid black histogram represents the fit result, (a) shows the AE 
projection for J-^jqq > 0.5 and r > 0.5. The blue hatched histogram shows the peaking background, 
the dashed red histogram shows the continuum contribution and the dash-dotted green histogram 
shows the total background, (b) shows the T^j^ projection for \AE\ < 0.01 GeV and r > 0.5. 
The dashed red histogram shows the continuum background contribution, (c) and (d) show the 
ni3 n and projections, respectively, for T^i^ > 0.5, \AE\ < 0.01 GeV and r > 0.5. The dashed 
blue histogram shows the af contribution, the red hatched histogram shows the contribution 
and the dash-dotted green histogram shows the total background. 
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T ai = 381 ± 43 (stat) MeV. From this fit, the relative contributions of each component are 
0.8% signal, 95.6% continuum, 3.3% BB background and 0.3% peaking background. 

Our measurement of the B° — > afir T branching fraction is lower than that measured by 
the BaBar Collaboration [l7[ by 1.9a though still in marginal agreement. We are also in 



agreement with predictions made in the QCD factorization framework given in Refs. [111. 112 



but not Ref. 13j. The upper limit of the branching fraction of B° — > afn T is also determined 
to be 

B(B° 4(1320)tt t ) x £(4(1320) 7r ± 7r T 7r ± ) < 2.2 x 10~ 6 at 90% CL, (19) 
which is improved over the current world average 24|| by about two orders of magnitude. 



D. Systematic Uncertainties 

Systematic errors from various sources are considered and estimated with independent 
internal studies and cross-checks. These are summarized in Table |V] This includes the 
uncertainty on the number of produced BB events in the data sample. Contributions to 
the uncertainty in the selection efficiency due to particle identification and tracking are 
calculated by independent studies at Belle. The systematic uncertainty arising from the 
assumption that the af decays exclusively through the dominant p°7r ± intermediate state is 
accounted for by recalculating the detection efficiency with an exclusive fa^OO)^ MC. 

The uncertainty in the af shape is determined by varying the fixed mass within its world 



average uncertainty [24J. We account for a difference in the fraction of mis-reconstructed 
signal events between data and MC by varying this parameter by ±5% and repeating the 
fit. Variations in the parametric model shape due to limited statistics are accounted for 
by varying each parameter within their errors. The dominant contribution to this category 
comes from the uncertainties in the signal shape correction factors obtained from analysing 
a high-statistics control sample B° — > D~[K + 7i~7t~]tt + . Uncertainties in the non-parametric 
shapes are obtained by varying the contents of the histogram bins within ±ler. The sys- 
tematic uncertainty due to fixing the peaking background yields is estimated by varying the 



branching fraction by its world average error [24j and repeating the fit. For modes where 
only an upper limit is known, the variation is taken as half of the upper limit. The fit bias 
is determined from pseudo-experiments by searching for a difference between the generated 
and fitted physics parameters. As channels containing an f are ignored in the nominal 
model, we account for a possible effect on the signal yield by embedding such events into 
these pseudo-experiments and determining further bias on the a\ir and ct27r branching frac- 
tions. Finally, the uncertainty from neglecting interference between a\ and a^ was estimated 
by constructing a 4-body amplitude including detector effects and generating three relative 
interference configurations between a\ and g^: maximum constructive interference, no inter- 
ference and maximal destructive interference. The largest deviation from the sample with 
no interference when fitting all with the nominal model gives the systematic uncertainty 
from interference. 
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TABLE V: Systematic uncertainties of branching fractions. 



Category 


5B(B° -> afirT) (%) 


6B(B° -> afirT) (%) 


N(BB) 


1.4 


1.4 


Tracking 


2.1 


2.1 


Particle identification 


4.5 


4.5 


/o(600) 


3.0 


N/A 


cti shape 


5.1 


19.4 


Mis-Reconstruction fraction 


2.1 


3.3 


Model shape 


4.2 


14.0 


Histogram shape 


4.8 


4.6 


Peaking background 


5.1 


11.4 


Fit bias 


0.4 


4.1 


Interference 


4.1 


4.1 


Total 


12.2 


28.1 



IV. TIME-DEPENDENT MEASUREMENT 



A. Event Selection 

In addition to the event selection criteria for the branching fraction measurement, events 
are selected for the time- dependent CP violation measurement if they satisfy |Ai?| < 
0.04 GeV. This requirement retains 97% of signal and 90% of peaking backgrounds, while 
rejecting 60% of the continuum background and 63% of the BB background. After this 
selection, the relative contributions of each component are 1.9% signal, 94.2% continuum, 
3.1% BB background and 0.8% peaking background. 



B. Event Model 



The signal PDF is given by 



V l £(At,q,c) = {1 + cAcp) 



-|A*I/Tb<> 

8t b o 



1 - qAw 1 ' 3 + q(l - 2w l ' s ) x 



(S C p + cAS) sin Am d At - (C C p + cAC) cos Am d At 



®R s B o BO (At), 



(20) 



which accounts for CP dilution from the probability of incorrect flavor tagging w l,s and the 
wrong tag difference Aw i,s between B° and B°, both of which are determined from flavor 



specific control samples [27J. This PDF is convolved with the At resolution function for 



neutral B particles R S B0 ^ , described in Ref. [25 



The reconstructed vertex position of mis-reconstructed events is dominated by the high 
momentum bachelor pion from B° — > afir T . This pion is rarely mis-reconstructed, and the 
effect of borrowing a lower momentum track from the tag-side results in a slightly smaller 
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lifetime. Thus, mis-reconstructed events are modelled using the truth model PDF with an 
effective lifetime and share CP parameters with the truth model. The effective lifetime is 
determined from mis-reconstructed signal MC events and the effect of this choice is accounted 
for in the systematic uncertainties. 

The continuum shape is determined from data taken below the T(4S) resonance with the 
model 

®i&(At). (21) 

This model contains a lifetime and prompt component to account for the charm and charm- 
less contributions, respectively, and is convolved with a sum of two Gaussians 

Rqq(At) = (1 — /tail)^(At; Mmean) ^main ") + /tail^(Ai; Mmeani ^main ^tail ) 7 (22) 

which uses the event- dependent At error constructed from the vertex resolution a = 
( \J °Rcc + ^Tag) / /^7 C as a scale factor. We also account for an asymmetry in the product 
qc, with the parameter AC q g, which is due to the jet-like topology of continuum. As a high 
momentum 7t + (/T~) in B^ cc is correlated with a high momentum 7r~(/T + ) on the tag-side, q 
and c will more often have the same sign. 

The BB shape is determined from a large MC sample and is divided into generic charmed 
and charmless, each further divided into neutral and charged BB decays, using a lifetime 
model 

V l - s (At,q,c) = l + C f tBB e — ®R BB (At). (23) 

4 ZT BB 

where R BB is relevant At resolution function for either neutral or charged B particles. 
Since reconstructed BB events may borrow a particle from the tag side, the average At 
lifetime tends to be smaller and is taken into account with the effective lifetime, t bb . Like 
continuum, BB events also exhibit a qc asymmetry; however, it is found to be more complex 
and is modelled with a 1st order polynomial in At 

AC BB = Po + Pl \At\. (24) 

The shapes of the peaking backgrounds are determined from individually generated MC 
events. As these backgrounds may exhibit CP violation, we use a model similar to the truth 
model, but with an effective lifetime. We also fix all time-dependent parameters to null with 
the effects of this choice reflected in the systematic uncertainties. 

To account for the broad underlying At events not yet described by either signal or 
background PDFs, a broad Gaussian outlier PDF is introduced 

V l 6 s ut (At,q,c) = ^G(At;0,a s Out ) 

The total likelihood for 83799 B° — > afn T candidates in the fit region becomes 

C = II II E -'T^'--- m ^ KJP'r^l'. q\ d) (26) 

l,S 1 = 1 j 
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P^(At,q,c) = 



1 + qcAC, 



[l-fs 



2t, 



fs S(At - /4) 



<-!<-! 



(25) 



TABLE VI: Statistical correlation matrix for the fit result. 



A C p Ccp AC S CP AS 



A C p 


1 








Ccp 


-0.20 


1 






AC 


+0.01 


+0.03 


1 




•Sep 


+0.02 


-0.01 


-0.03 


1 


AS 


+0.09 


-0.05 


-0.01 


+0.02 



where /' s is the event-dependent probability of component j, in flavor-tag bin /, with de- 
tector configuration s 



h ^ bb/gg ,m 3n , n 37T ) - N sfhs v i,s (AF i T i m i V i y W> 



bb/qq' '"37T! ' l 3ttJ ~ ^ ATS f l,s^.l,s 



constructed from the branching fraction measurement. Only the 5 time-dependent coeffi- 
cients of the truth model are free in the fit to data. 

As a consistency check, we perform a fit to data to measure the B° lifetime while fixing 
the 5 parameters of the truth model to zero. We obtain r^o = 1.389 ± 0.085 ps, which is in 
agreement with the current world average r^o = 1.519 ± 0.007 ps 24 



C. Fit Result 

We perform a fit to the data and obtain the CP violating parameters 

A CP = -0.06 ± 0.05 (stat) ± 0.07 (syst), 

C CP = -0.01 ± 0.11 (stat) ± 0.09 (syst), (28) 
S CP = -0.51 ± 0.14 (stat) ± 0.08 (syst), 

and the CP conserving parameters 

AC = +0.54 ± 0.11 (stat) ± 0.07 (syst), 

AS = -0.09 ± 0.14 (stat) ± 0.06 (syst). ^ ' 

The background subtracted fit results are shown in Fig. HJ where the data points are the sig- 
nal yields obtained by repeating the branching fraction fits in At, q or At, qc bins, accounting 
for the selection criteria on AE. Ours are the most precise measurements of these param- 



eters to date and are in agreement with those obtained by the BaBar Collaboration [18 



We are also in very good agreement with the theoretical predictions made in Ref. [111 ]. The 
statistical correlation coefficients between the obtained parameters are given in Table I VI I 

A MC technique is employed to obtain ifi^ using Eq. H] in order to take correlations be- 
tween the fitted parameters into account. We generate multiple vectors (Ccp, AC, Sep, AS) 
based on a correlated mult i- dimensional Gaussian constructed from the fit result and solve 
for the four solutions of 0| ff each time. We take the central values and their uncertaintes 
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FIG. 4: (color online) Background subtracted time-dependent fit results for B° — > a 1 7r T . (a) and 

(b) show the At distributions for the -E?Tag flavor q and the product of the fl avor an d ai charge 
qc, respectively. The left plots are useful for visualizing the effect of flavor-dependent CP violation 
while the plots on the right show the effects of the CP conserving parameters. The solid blue and 
dashed red curves represent the At distributions for positive and negative quantity, respectively. 

(c) and (d) show the asymmetry of the plots immediately above them, (iV^ 1 * — A r ^ lt )/(A^^ lt + A r ^ lt ), 
where Nj} 1 (iVE 1 *) is the measured signal yield of positive (negative) quantities in bins of At. 

from the resulting distributions of (pf 1 and obtain the four solutions 

0f = (-17.3 ± 6.6 (stat) ± 4.8 (syst)) , 
<f>f = (41.6 ± 6.2 (stat) ± 3.4 (syst)) , 
<j>f = (48.4 ±6.2 (stat) ± 3.4 (syst)) , 

0f = (107.3 ± 6.6 (stat) ± 4.8 (syst))°. (30) 
Using a similar technique, we obtain the direct CP violation parameters given in Eq. [5] 

A + _ = +0.07 ± 0.08 (stat) ± 0.10 (syst), 

A_ + = -0.04 ± 0.26 (stat) ± 0.19 (syst), (31) 

where the statistical correlation coefficient between these two parameters is 0.61. 

We compose Gaussian distributions from the four solutions for 0| ff and construct a two- 
sided p- value plot for c^| ff as shown in Fig. [5j We also perform a likelihood scan to estimate 
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FIG. 5: Difference 1-CL plotted for a range of </>| values as shown by the solid curve. The p-value 
is taken as the maximum p- value that can be calculated from all four solutions. The dashed line 
indicates CL = 90%. 




S CP AC 

FIG. 6: (color online) Likelihood scan of Sep (a) and AC (b). The dotted blue curve shows the 
statistical likelihood while the solid black curve includes the systematic uncertainty. 

the significance of Sep and AC as shown in Fig. EJ The significance of mixing-induced CP 
violation is found to be 3.1a including systematic uncertainties while the rate where the af 
does not contain the spectator quark is found to dominate the rate where it does at the 4.1a 
level. 

D. Systematic Uncertainties 

The systematic uncertainties for the time- dependent parameters are summarized in Ta- 
ble I VIII These are estimated from various sources including uncertainties in the IP profile, 
charged track selection based on track helix errors, helix parameter corrections, At and ver- 
tex goodness-of-fit selection, Az bias and SVD misalignment. The fixed physics parameters 
7~£0 and Arrid, At resolution function and data model shape parameters including background 
effective lifetimes and qc asymmetries, as well as the flavor tagging performance parameters 
w and Aw, are varied by ±la. We generate MC pseudo-experiments and perform an en- 
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semble test to obtain systematic biases from interference on the tag-side arising between the 
CKM favored bd — > (cud)d and doubly-CKM-suppressed bd — > {ucd)d amplitudes in final 
states used for flavor tagging [30j . These sources should not affect Acp, as this parameter 
represents a time and flavor-integrated asymmetry. 

The remaining systematic categories affect all time- dependent parameters. The parame- 
ters and non-parametric shapes describing signal probability are varied in the same way as 
for the branching fraction measurement. The CP violation parameters of mis-reconstructed 
events are assumed to be the same as signal. To account for this, a sample of GEANT MC is 
produced with the nominal fit result. The systematic error is taken as the difference between 
the fit result to the correctly reconstructed subsample and a fit to the whole sample sharing 
the CP parameters between the signal and mis-reconstructed components. The fit bias is 
determined from an ensemble test by searching for a difference between the generated and 
fitted physics parameters. 

Possible CP violation in the background is the dominant systematic uncertainty. We 
assume that the neutral BB background possesses a 20% CP violating component and that 
peaking backgrounds display a 50% CP violating effect. For the neutral BB component, 
the uncertainty is determined by fixing the CP parameters and refitting the data. For the 
peaking backgrounds, special GEANT MC samples are produced with the nominal fit result 
for signal and CP violation generated in the peaking background. We perform a fit with 
the peaking background CP parameters fixed to null asymmetry and compare this with a 
fit where they are fixed to the generated values. 

To estimate the effects of interference on the signal probability, we employ the same 4- 
body amplitude generator described in the branching fraction measurement. MC samples 
are generated with the nominal fit result with the three phase configurations: maximum 
constructive interference, no interference and maximal destructive interference. Firstly, sig- 
nal probability is predetermined for all three samples by fitting for the branching fractions 
of the a\Ti and a27r components assuming no interference in the fit model. Finally, a time- 
dependent fit is performed to the three phase configurations and the systematic error taken 
as the maximum deviation from the sample with no generated interference. 



V. CONCLUSION 

We have presented a measurement of the product branching fraction and time-dependent 
parameters in B° — > af(12Q0)iT T decays, which are in agreement with measurements per- 
formed by the BaBar Collaboration 17, 3]. We obtain the product branching fraction 



B(B° a^(1260)7r T ) x B(af (1260) ttVV 1 ) = (11.1 ± 1.0 (stat) ± 1.4 (syst)) x 10~ 6 , 

(32) 

and an upper limit on the product branching fraction for a possible decay with the same 
final state 

B(B° 4(1320)tt t ) x £(af(1320) ttVtt*) < 2.2 x 10~ 6 at 90% CL. (33) 

This upper limit is an improvement over the current world's most restrictive limit by about 
two orders of magnitude. In a time-dependent measurement to extract CP asymmetries, we 
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TABLE VII: Systematic uncertainties of time-dependent parameters. 


Categories 


related to At 


reconstruction and flavor tag 


ging are not applicable to the time and flavor integrated Acp- 


Category 


SA C p (ltr 2 ) 


SC C P (10- 2 ) 


SAC (10~ 2 ) 


SS C p (10- 2 


) SAS (10- 2 ) 


IP profile 


N/A 


0.2 


0.2 


1.0 


1.0 


-^Tag track selection 


N/A 


1.2 


0.4 


0.8 


1.1 


Track helix errors 


N/A 


0.0 


0.0 


0.0 


0.0 


At selection 


N/A 


0.1 


0.0 


0.1 


0.0 


Vertex quality selection 


N/A 


0.3 


0.9 


0.3 


0.2 


Az bias 


N/A 


0.5 


0.5 


0.4 


0.4 


Misalignment 


N/A 


0.4 


0.4 


0.2 


0.2 


t b o and Am c i 


N/A 


0.3 


0.3 


0.2 


0.2 


At resolution function 


N/A 


1.3 


0.9 


2.8 


1.7 


Flavor tagging 


N/A 


0.3 


0.2 


0.2 


0.1 


Model shape 


N/A 


0.3 


2.9 


0.6 


0.5 


Tag-side interference 


N/A 


3.6 


0.2 


0.5 


0.4 


Signal probability 


0.5 


0.4 


1.9 


2.0 


0.8 


Mis- Reconstruction 


0.2 


0.1 


0.7 


0.3 


0.3 


Fit bias 


0.4 


1.6 


0.3 


1.0 


0.1 


Background CP violation 


6.6 


7.6 


5.7 


6.8 


5.1 


Interference 


0.8 


1.1 


2.0 


0.2 


1.0 


Total 


6.6 


8.9 


7.2 


7.9 


5.8 



obtain the CP violation parameters 

A CP = -0.06 ± 0.05 (stat) ± 0.07 (syst), 

C CP = -0.01 ± 0.11 (stat) ± 0.09 (syst), (34) 
S CP = -0.51 ± 0.14 (stat) ± 0.08 (syst), 

representing time and flavor integrated direct, flavor-dependent direct and mixing- induced 
CP violation, respectively. Simultaneously, we also extract the CP conserving parameters 

AC = +0.54 ± 0.11 (stat) ± 0.07 (syst), 

AS = -0.09 ± 0.14 (stat) ± 0.06 (syst), ^ > 

which, respectively, describe a rate difference and strong phase difference between the decay 
channels where the af does not contain the spectator quark and those where it does. We 
find first evidence of mixing- induced CP violation in B° — > a^(1260)7r T decays with 3.1a 
significance and the rate where the af does not contain the spectator quark is found to 
dominate the rate where it does at the 4. la level. However, there is no evidence for either 
time and flavor integrated direct CP violation or flavor-dependent direct CP violation. Our 
results are in good agreement with theoretical predictions given within the QCD factorization 
framework (ill. 121 ] and may be used in an either an isospin analysis [l5|] or SU (3) flavor 



symmetry [14J to extract 
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